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K0  modified  Bessel  function 

J^(X)  integrated  kernel  function,  see  Equation  (3.8) 
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I.  INTRODUCTION 


In  Reference  1  the  viscous  theory  of  unsteady  lift  on  two-dimensional 
airfoils  of  arbitrary  thickness  distribution  is  developed  in  detail. 
EAtensive  numerical  examples  are  given  for  steady  and  unsteady  flow  that 
compare  favorably  with  experimental  results.  Both  the  d'Alembert  paradox  and 
the  fatness  paradox  (see  Reference  2)  are  resolved,  and  most  important,  the 
role  of  viscosity  in  establishing  a  unique  value  of  the  circulatory  part  of 
the  lift  is  completely  understood  within  the  framework  of  this  theory.  The 
concept  of  viscous  thin  airfoil  theory  is  directly  applicable  to  other  lifting 
bodies,  e.g.,  three-dimensional  wings  and  bodies  of  revolution.  In  a  recent 
report  (Reference  3),  completed  under  contract  to  NASA,  the  preliminary 
application  of  the  viscous  theory  to  a  rectangular  wing  is  presented.  The 
numerical  results  are  encouraging  and  reinforce  the  basic  principles  of  the 
two-dimensional  theory. 

The  present  report  is  directed  towards  the  application  of  the  viscous 
theory  to  bodies  of  revolution.  The  complete  unsteady  potential  theory  is 
presented  and  subsequently  corrected  for  the  effect  of  viscosity.  The  theory 
is  applied  to  the  simple  case  of  a  ring-wing.  Both  numerical  results  and 
asymptotic  results  are  presented.  For  very  slender  ring-wings  we  recover  the 
well  known  results  of  slender  body  theory.  For  very  large  diameter  to  length 
ratios  we  obtain  a  simple  asymptotic  result  that  is  in  exact  agreement  with  a 
theoretical  result  of  Ribner  (Reference  4).  Fhe  theoretical  result  is  also  in 
approximate  agreement  with  experimental  results  of  Flatau  (Reference  5). 
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II.  VISCOUS  THEORY  OF  LIFT  FOR  BODIES  OF  REVOLUTION 
IN  INCOMPRESSIBLE  FLOW 
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A.  Preliminary  Remarks 

Consider  a  body  of  revolution  in  a  high  Reynolds  number  (Re  >  10®) 
incompressible  flow  (see  Figure  2.1).  In  the  absence  of  any  disturbance  that 
would  produce  a  transverse  force  (lift),  the  body  axis  is  aligned  with  the 
free  stream  and  the  mean  steady  flow  is  characterized  by  the  velocity  field  vQ 
that  we  assume  to  be  given.  Suppose  that  the  body  with  its  boundary  layer  and 
wake  are  enclosed  by  a  cylindrical  control  surface  S.  In  the  domain  D 
exterior  to  S  the  flow  is  free  of  vorticity  and  thus  has  a  velocity  potential. 
Forward  of  the  plane  P  in  Figure  2.1,  the  distance  between  the  surface  S  and 
the  actual  body  surface  will  be  of  the  order  of  the  boundary  layer  thickness. 
The  plane  P  can  be  imagined  to  lie  a  few  percent  of  the  body  length  forward  of 
the  aft  end.  Immediately  downstream  of  the  body,  the  surface  S  is  faired 
smoothly  into  a  cylinder  of  constant  radius  that  bounds  the  wake.  As  the 
Reynolds  number  is  increased  beyond  all  bound  the  radius  of  the  wake  cylinder 
will  tend  to  a  constant. 

Forward  of  the  plane  P,  the  normal  momentum  and  pressure  at  a  point  on  S 
will  differ  from  the  corresponding  values  on  the  actual  body  surface  by  terms 
that  are  of  the  order  6  /L  where  *>  is  the  turbulent  boundary  layer 
displacement  thickness  and  L  is  the  local  radius  of  curvature  or  the  body 
length,  whichever  is  smaller.  Terms  of  this  order  will  be  neglected  in  the 
present  theory.  Downstream  of  the  body  the  transverse  load  on  S  must  tend  to 
zero. 


To  derive  the  viscous  theory,  the  first  step  is  t^  formulate  the 
potential  flow  problem  in  D  and  then  recast  it  in  the  form  of  an  integral 
equation  for  the  transverse  load  function  (to  be  defined),  given  that  the 
normal  momentum  (upwash)  on  S  is  specified.  The  local  singularity  in  the 
load-upwash  kernel  function  is  then  modified  (reduced  in  strength)  to  account 
for  the  direct  effect  of  viscosity.  The  principle  for  introducing  the  local 
effect  of  viscosity  is  based  on  our  previous  work  with  viscous  wing  theory 
(Reference  1). 
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B.  Potential  Flow  Problem 


Let  $  denote  the  perturbation  potential  due  to  an  unsteady  simple 
harmonic  transverse  momentum  disturbance  on  S.  The  boundary  value  problem  for 
4  is  the  following: 


V2<j>  =  0  in  D 


(2.1) 


M  .  /fr)sine  on  S:  r  =  R(x) 


(2.2) 


grad  $  ~  o  at  »  (except  near  the  wake) 
and  for  the  perturbation  pressure 


(2.3) 


P  = 


(2.4) 


where 


(2.5) 


and  u  is  the  radian  frequency  of  the  specified  transverse  momentum 
disturbance.  Also  note  that 


-  R‘ 


JL 

a* 


JU  JL 
ax  ar 


(2.6) 


where 


n  =  -  R' (x)r  +  t 


(2.7) 


is  normal  to  the  surface  S.  The  unit  normal  on  S  is  given  by 


|fl|  \/l  +  R'2 


(2.8) 


We  use  Cartesian  and  cylindrical  coordinates  interchangeably  in  the  subsequent 
development. 


Before  we  proceed  with  the  representation  of  <|>,  we  first  consider  the 
inverse  of  the  relation  (2.4).  Note  that 


.►  5  d 

V°  *  ^  =  q°  d? 


(2.9) 


where 


ds  =  s'dx  =  +  R*  2  dx 


%  =  lvo 


the  speed  of  the  nonlifting  flow.  Thus,  we  can  write 


84  ,  1»  A  _ 
ds  90  *  " 


(2.1U) 


(2.11) 


(2.12) 


The  solution  of  this  equation  for  ♦  is  straightforward.  The  final  result  is 
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34%?S4 


♦(x,S)  =  v.  f  X(5,0)e  -i“T(x,5)  d!r 


(2.13) 


where 


f»/l+Rl2 

p®q0v» 


(2.14) 


T(x,5)  = 


fl+R*2 


(2.15) 


We  refer  to  X  as  the  generalized  load  function  since  it  defines  completely  the 
local  lift  force  on  the  surface  S.  The  two  point  function  T(x,£)  is  the  time 
required  for  a  vortical  disturbance  produced  by  a  harmonic  load  at  station  5 
to  convect  along  the  surface  at  speed  qQ  to  station  x.  It  is  the  most 
important  unsteady  effect  of  body  shape  on  the  aerodynamic  loads. 

C.  The  Load-Upwash  Integral  Equation 
Consider  the  Green's  function 


(2.16) 


that  satisfies  the  Poisson  equation 


V2G  =  -  4 it  6(5 ?-$) 


(2.17) 


where  6(5?)  is  the  Dirac  delta  function.  Form  the  identity 


GV2*  _  $y2g  =  4n  4>  6(5?-$) 


(2.18) 
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and  integrate  with  respect  to  coordinates  %  over  the  domain  D.  After  the 
application  of  the  divergence  theorem,  the  final  result  can  be  expressed  in 
the  following  form: 


♦  (*) 


_1_ 

4n 


(2.19) 


where 


dA  =  R(5)d£  d0' 


(2.20) 


and 


a 

aN ' 


ap 


(2.21) 


The  next  step  is  to  eliminate  the  polar  angle  9  from  the  problem.  It 
follows  from  the  boundary  condition  (2.2)  that  all  dependent  variables  are 
proportional  to  sin  9.  Thus,  we  write 

4>  «  5*(x,r)sin  0  (2.22) 

P  =  P(x,r)sin  9  (2.23) 

and  also  for  the  load  function 

X  =  .§tx,r)sin  9  (2.24) 

Substitute  Equations  (2.22)  and  (2,2)  into  (2.19)  to  obtain 
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where 


*(x,r) 


_ cosjje _ 

\/x2+r2+p2-2rp  cos  6 


4 

Vx2+(r+p)2 


K(k)  +  4  E(k) 

kz 


with 


k2  _  4rp 

x2+(r+p)2 


E(k) 


-*/2  - - 

=  J  yl-k2sin* 


e  dQ 


K(k)  = 


/**/2 

Jo 


de 

\/l-k2sinz0 


The  functions  E(k)  and  K(k)  are  standard  elliptic  integrals  (see 
Now  replace  ®  in  the  first  term  of  (2.25)  with  the  integral 
function;  i.e.. 


(2.25) 


(2.26) 


(2.27) 


(2.28) 


(2.29) 


(2.30) 


Reference  6). 
over  the  load 
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*(5,R(5))  •  v 


(2.31) 


-1»T(W)  d5. 


The  result  is 


®(x,r) 


V=o 


/. 


e  -1uTU'.5) 


atT 


R(5* )d5' 


4n 


Finally,  we  compute  the  normal  derivative  of  the  last  result 
surface  S.  The  last  term  has  a  Cauchy  principle  value  so  that 


Lim 

r-*-R(x) 


and  the  slash  through  the  integral  indicates  the  usual  limit  operation 


d£  =  Lim 
e+0 


d5 


The  integral  equation  forj^can  thus  be  written  in  the  following  form: 


■■  -f  ]S?U)%{x,l)K  =>^x)/V„ 

JxN 


(2.32) 

on  the 

(2.33) 

i.e., 

(2.34) 

(2.35) 


where 


Km 


-  KA  *  j:  f 

J X, 


(2.36) 


*rLJ?x)  /  e',“T(5,.5)  ^pr*(fW  (2.37) 


The  upper  limit  of  integration  in  Equation  (2.35)  is  Xy  since  the  transverse 
load  must  vanish  on  the  wake  cylinder  (to  the  order  we  are  working).  The 
generalized  upwash  /^(x)  can  be  separated  into  two  parts;  i.e.. 


Km  -K'm  *£1m 


(2.38) 


where 


(2.39) 


"i, .  i  r  & 

i  (X)  =2^  /  5N 
*'XT 


(2.40) 


The  first  part,  /^%  is  completely  determined  in  terms  of  the  known  upwash  on 
the  surface  of  the  body.  The  second  part,  can  only  be  determined  after 
)! is  calculated  on  the  wake  cylinder.  The  actual  numerical  process  must  be 
carried  out  iteratively.  First  assume  >^00  to  be  zero  and  >olve  Equation 
(2.35)  for  the  load  function.  Then  ca’culate  on  the  wake  cylinder  by  using 
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Equation  (2.35)  with  x  >  Xy.  Then  calculate  ftQ-  with  Equation  (2.40)  and 
correct  the  generalized  upwash  with  Equation  (2.38).  If  the  wake  correction 
is  shown  to  be  small,  the  process  can  be  terminated.  Otherwise,  it  will  be 
necessary  to  recalculate  a  new  load  function. 

We  emphasize  that  the  foregoing  problem  is  not  well  posed  as  yet.  The 
reason  is  that  the  kernel  function  (2.37)  has  been  derived  with  potential 
theory  and  has  a  Cauchy  singularity  (see  below  Equation  (2.46)).  The  integral 
equation  (2.35)  has  an  eigensolution.  The  purpose  of  the  viscous  theory  is  to 
correct  the  singularity  in  the  kernel  such  that  Equation  (2.35)  has  a  unique 
solution.  We  turn  now  to  an  explicit  evaluation  of  the  kernel  for  the  case  of 
a  ring-wing.  The  reason  is  to  illustrate  the  basic  structure  of  the  kernel 
and,  in  particular,  the  singularity.  Also  we  show  how  conventional  slender 
body  theory  "tampers"  with  the  singularity  in  such  a  way  that  a  unique 
solution  is  obtained  in  terms  of  an  effective  base  area. 

D.  Properties  of  the  Kernel  Function  for  a  Ring-Wing  in  Steady  Flow 

Consider  the  ring-wing  shown  in  Figure  2.2.  For  to  =  0  the  kernel  (2.37) 
can  be  expressed  in  the  following  form: 

=  a^fx-5)  (2.41) 

where 


jn*) 


(2.42) 


The  detailed  evaluation  of  the  ring-wing  kernel  is  given  in  the  Appendix.  The 
final  result  is 


t  f +  [®r  ft ' k)2  ‘  ?  (1-kZ)3/2  K(k)] S9n  x 


(2.43) 


where 


k  =  TfN  (2*44) 

^x2+4a2 

and  E(k),  K(k)  are  the  elliptic  integrals  given  by  Equations  (2.29)  and 
(2.30).  Also  we  have  the  following  two  asymptotic  results  for  large  and  smal t 
x,  respectively: 

For  x  >  2a 

JT(x)  ^  -  ~  H(x)  t?.45) 

a* 

For  x  «  2a 

M*)  =  -  f  *  7  (2.46) 

Numerical  calcul  itions  of  the  load  distribution  on  a  ring-wing  may  be 
found  in  Section  Hi.  Here  we  want  to  synthesize  the  viscous  correction  to 
the  kernel  that  is  necessary  to  obtain  a  unique  solution  of  the  integral 
equation.  Also  we  can  compare  the  viscous  and  inviscid  kernels  with  the 
asymptotic  kernel  (2.45)  that  yields  the  results  of  slender  body  theory. 

E.  The  Viscous  Slender  Body  Kernel  Function 

The  asymptotic  behavior  of  the  kernel  function  for  small  x  (i.e., 
x  «  2a,  see  Equation  (2.46))  dictates  a  fundamental  mathematical  property  of 
the  integral  equation  (2.35).  That  is,  the  equation  has  an  eigensolution 
whose  magnitude  is  proportional  to  the  lift  on  the  body.  The  problem  is 
completely  analogous  to  the  two-dimensional  wing  theory  discussed  in  Reference 
1.  There  we  have  shown  rigorously  that  the  correct  viscous  modification  of 
the  kernel  is  to  convert  the  Cauchy  singularity  into  a  weaker  logarithmic 
singularity  on  the  Stokes  scale.  For  the  ring-wing,  the  correct  modification 
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of  the  asymptotic  result  (2.46)  is 


M*)  S  '  7  +  eV"X/2v  K0(vjx|/2v)j 


(2.47) 


where  KQ  is  the  modified  Bessel  function.  For  x  «  2a  the  viscous  ring-wing 
kernel  behaves  like 


=  7  tn  |  x  j  ,  x  «  2a 


(2.48) 


The  argument  for  correcting  the  ring-wing  kernel  for  viscosity  is  based 
entirely  on  our  previous  work  on  the  two-dimensional  wing.  So  long  as  the 
effect  of  viscosity  is  confined  to  a  circumferential  region  that  is  small 
compared  to  the  smallest  geometric  dimensions,  the  ring-wing  behaves  locally 
like  a  two-dimensional  wing.  If  we  introduce  the  following  dimensionless 
parameters. 


(2.49) 


(2.50) 


then  a  necessary  requirement  for  the  validity  of  Equation  (2.47)  is  that 


(2.51) 


This  result  says  that  the  Stokes  viscous  length  scale  must  be  much  smaller 
than  the  smallest  geometric  length  scale. 


I 


U 


Another  remarkable  property  of  the  kernel  function  is  that  the 
singularity  near  the  load  is  not  changed  when  the  flow  become:,  unsteady.  Also 
the  complete  unsteady  potential  theory  kernel  for  an  arbitrary  body  of 
revolution  (see  Equation  (2.37))  ruts  a  local  Cauchy  singularity  like  a 
ring-wing  of  radius  R(x).  Thus,  we  can  write  down  the  complete  viscous 
slender  body  kernel  without  further  ado: 


*  rL4?x> 


f. 


,-iMT(5‘  ,£) 


3N5N* 


Hl'W' 


(2.52) 


2  v.(x-5)/2v 


R(x) 


K0(v»|x-C|/2v) 


The  complete  problem  of  calculating  tne  lift  (steady  or  unsteady)  on  a  body  of 
revolution  is  to  solve  the  integral  equation  (2.35)  with  the  viscous  kernel 
(2.52).  The  solution  is  unique  because  the  kernel  has  the  correct  logarithmic 
singularity.  Specific  calculations  for  a  ring-wing  are  given  in  Section  III. 

F.  The  Kernel  of  Slender  Body  Theory 

We  assume  that  the  asymptotic  result  (2.45)  is  valid  for  all  x  and 
furthermore,  extend  it  to  a  slender  body  of  revolution  by  replacing  a  with 
R(x).  Then  the  integral  equation  (2.35)  can  be  written  in  the  form 


2n 

R2(x) 


-H(x-5)  «-*,(«) 


(2.53) 


where 


JTt)  - 


(2.54) 
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I 


I 
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I 


is  local  the  force  per  unit  length  on  the  body.  The  solution  of  Equation 
(2.53)  is 


J*Tx) 


(2.55) 


For  a  body  at  constant  angle  of  attack,  we  have 


(2.56) 


and 


I 


(2.57) 


* 


where 


S  =  nR2 


(2.58) 


*  is  the  local  cross-sectional  area  of  the.  body.  We  recognize  the  last  result 

as  the  familiar  formula  of  classical  slender  body  theory;  i.e.,  the  local 
lift  force  on  a  body  of  revolution  at  constant  angle  of  attack  is  proportional 
to  the  derivative  of  the  cross-sectional  area.  The  slender  body  formula 

»  yields  a  reasonably  good  approximation  of  the  sectional  lift  on  the  forward 

portion  of  the  body.  Unfortunately,  for  closed  bodies  of  revolution,  it  leads 
to  the  embarrassing  result  that  the  total  lift  is  zero.  In  application  of  the 
theory  it  has  become  cumstomary  to  introduce  an  "effective  base  area"  and  to 
correlate  measured  lift  data  in  terms  of  this  artificial  parameter.  With  the 
foregoing  general  theory  we  can  offer  a  plausible  explanation  of  this  apparent 
uniqueness  that  is  brought  about  by  introducing  the  familiar  slender  body 

approximation  into  the  potential  flow  theory.  We  use  the  term  apparent 

uniqueness  because  the  unknown  lift  is  simply  traded  off  for  an  equally 

unknown  base  area. 
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The  potential  flow  kernel,  the  viscous  kernel  and  the  slender  body 
aporoximation  thereof  are  plotted  in  Figure  2.3.  From  the  point  of  view  of 
the  theory  of  integral  equations,  the  slender  body  approximation  represents  a 
rather  crude  elimination  of  the  fundamental  Cauchy  singularity  in  the  exact 
potential  flow  theory.  Indeed,  this  step  eliminates  the  possibility  of  an 
eigensolution  and  thus  yields  a  unique  value  for  the  local  normal  force.  The 
approximation  is  not  too  bad  but  does  not  lead  to  a  unique  value  for  tne  total 
lift.  The  viscous  modification  of  the  Cauchy  singularity  is  the  correct  way 
to  proceed.  A  unique  value  of  the  total  lift  can  be  calculated  without  tne 
device  of  an  artificial  base  area. 


) 


i 


> 


* 


* 
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III.  NUMERICAL  RESULTS  FOR  A  RING-WING  AIRFOIL  IN  STEADY  FLOW 

To  calculate  the  load  distribution  on  a  ring-wing,  it  is  convenient  to 
introduce  dimensionless  variables  as  follows: 

X  =  ~  (3.1) 


4a  _  diameter 
A  half  length 


(3.2) 


k  =  t/  x2  +  T2 


(3.3) 


vjt 

o  =  — —  Reynolds  Number  referred  to  A/4 
4v 


(3.4) 


If  we  choose  the  origin  of  coordinates  to  be  on  the  axis  at  the  midpoint  on 
the  body  axis,  then  the  normalized  integral  equation  can  be  written  in  the 
following  form: 


_L 

Zn 


J& 


^V^t"X-Y)dY  *  «(X) 


(3.5) 


where 


a(X)  =  -  z0'  local  angle  of  attack  of  the  body  axis  (3.6) 


jn«)  -  (3.7) 
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and 


MM  =  T  -  kW)  +  e0X  K0(a(x|)  (3.8) 

8  f1* 

+  -%  /  (|x|  -2s)T(k')ds 

xZ  Jo 


k’  =  T/  Js2+t2  (3.10) 

We  refer  toJ^(x)  as  the  integrated  kernel. 

The  overall  normal  force  and  moment  coefficients  have  the  following 
integral  representations: 


■/; 


J^U)dX  = 


2)  H 


(3.11) 


(3.12) 


The  reference  area  is  chosen  to  be  one  half  the  surface  area  of  the  ring-wing. 
The  moment  is  measured  positive  nose  up  about  the  center  of  the  ring. 
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The  normalized  integral  equation  (3.5)  is  precisely  of  the  same  form  as 
the  equation  of  two-dimensional  viscous  thin  airfoil  theory  (see  Reference  1, 
Equation  (3.25)).  The  kernel  has  a  logarithmic  singularity  and  this  renders  a 
unique  solution  of  Equation  (3.5).  The  load  function  must  admit  square  root 
singularities  at  the  fore  and  aft  ends.  The  equation  was  reduced  to  a  finite 
spectral  form  with  an  assumed  polynomial  approximation  of  the  load;  i.e.. 


where  Tp  is  the  Chebyshev  polynomial  of  the  first  kind  of  degree  n  (see 
Reference  6).  The  numerical  algorithm  is  discussed  in  detail  in  Reference  1. 
The  integral  equation  was  solved  for  a  computational  Reynolds  number  of  1000 
with  20  Chebyshev  polynomials  to  approximate  the  load  function.  It  was 
determined  that  the  solution  was  adequately  converged  for  this  choice  of 
parameters. 

Typical  calculations  of  the  load  distribution  are  plotted  in  Figures  3.1, 
3.2  and  3.3  for  three  values  of  the  fineness  ratio  x.  We  note  from  these 
calculations  that  the  normal  force  coefficient  is  proportional  to  x  and  thus 
tends  to  zero  as  x  tends  to  zero.  This  result  can  also  be  deduced 
analytically  since  the  kernel  is  proportional  to  1/x  for  very  slender  bodies. 
The  limiting  form  is  in  agreement  with  the  result  of  slender  body  theory  that 
predicts  a  zero  lift  force  on  any  u.dy  of  constant  radius.  For  larger  values 
of  the  fineness  ratio  (see  Figure  3.1  with  x  =  1)  the  load  distribution  tends 
to  a  form  that  is  similar  to  the  lift  on  a  two-dimensional  airfoil.  This 
result  can  also  be  deduced  analytically  from  the  asymptotic  form  of  the  kernel 
(see  Equation  (2.46))  and  is  intuitively  expected  since  each  section  of  a  thin 
hoop  acts  like  a  two-dimensional  airfoil  section  at  a  local  angle  of  attack 
that  is  a  maximum  on  the  upper  and  lower  section  and  tends  to  zero  on  the 
sides. 

To  illustrate  the  last  point  we  note  that  the  normalized  ring-wing 
problem  can  be  expressed  in  the  following  form  for  x  »  l: 


I  f  » 


J  Sfr)  ^  j*i|X-Yl  +  e°(x-Y)  K0(o|X-Y|)J 


dY  =  a  -  ctj  (3.14) 


where 


^lY)dY  =  ^ 


(3.15) 


In  the  limit  a  ■*■  <»,  Equation  (3.13)  has  the  unique  solution. 


& (x)  ‘V==? 


(3.16) 


and  if  we  integrate  the  last  result  over  the  chord,  we  obtain 


r  _  tt 
N  ~  1+h/t 


(3.17) 


Ring-Wing 


1 

1+k/t 


2-D  Wing 


(3.18) 


The  last  result  is  precisely  the  same  form  as  the  lift  curve  slope  of  a 
three-dimensional  wing  with  the  typical  aspect  ratio  correction.  It  also 
agrees  exactly  with  a  result  derived  by  Ribner  (Reference  4)  with  the  more 
elementary  lifting  line  theory. 


We  compare  the  last  asymptotic  result  with  the  experimental  data  of 
Flatau  (Reference  5)  for  a  11.7%  thick  Clark  Y  airfoil.  Typical  section  data 
are  presented  in  Figure  3.4  for  i  *  4.9.  From  the  CL  versus  o  data  we 
estimate  the  lift  curve  slope  to  be  2.0/rad.  From  Equation  (3.17)  the 
theoretical  result  is  2.8/rad,  a  value  somewhat  larger  than  the  experimental 
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result.  The  discrepancy  is  probably  due  to  two  sources.  First,  the 
theoretical  formula  is  weakly  asymptotic  in  the  parameter  it/x.  For  example, 
if  the  next  term  in  the  asymptotic  expansion  of  the  denominator  of  Equation 
(3.17)  has  an  additive  term  of  order  (*/t)2  then  for  t  =  4.9  (the  Flatau 
value),  we  can  expect  a  25%  reduction  in  the  theoretical  value.  Also  the 
two-dimensional  section  lift,  coefficient  of  a  Clark  Y  airfoil  is  5.25/rad 
(Reference  7)  compared  to  the  flat  plate  value  of  2n.  This  represents  another 
16%  reduction  in  the  calculated  theoretical  value. 
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IV.  SUMMARY  AND  RECOMMENDATIONS 


The  viscous  theory  of  lift  on  bodies  of  revolution  has  been  derived  in 
detail.  The  theory  is  reduced  to  an  integral  equation  of  the  airfoil  type 
with  a  kernel  function  that  has  a  logarithmic  singularity.  Detailed 
calculations  of  the  load  distribution  and  several  asymptotic  results  are 
presented  for  a  ring-wing  airfoil  (i.e.,  the  simplest  body  of  revolution). 
The  calculated  load  distributions  are  similar  to  those  on  a  two-dimensional 
flat  plate  airfoil  for  large  diameter  to  length  ratios.  It  is  further  shown 
that  the  asymptotic  normal  force  coefficient  on  a  very  large  diameter 
ring-wing  is  in  agreement  with  a  theoretical  result  of  Ribner  (Reference  4). 
For  very  slender  ring-wings  the  normal  force  coefficient  tends  to  zero  in 
agreement  with  the  prediction  of  slender  body  theory. 

It  is  recommended  that  the  complete  steady  flow  theory  for  general  bodies 
of  revolution  be  programmed  and  numerical  results  obtained  for  the  load 
distribution.  The  theoretical  calculations  should  be  compared  to  known 
experimental  data. 
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APPENDIX 

Evaluation  of  the  Ring-Wing  Kernel  (Potential  Flow) 


Consider  the  following  kernel  function: 


vys.  I, 


(A.l) 


where 


cos  e' 


(A.2) 


R'  =  52  +  r2  +  p2  _  2 rp  cosQ' 


(A. 3) 


Note  that 


i  =  ^  5  .  p-r  cos8‘ _ 1 

R  5rd^  [_  R  r2+p2-2rp  cos8'J 


(A.4) 


Introduce  the  change  of  variable 


8*  =  it  -  6 


(A. 5) 


and  substitute  Equation  (A.4)  into  Equation  (A.l).  The  final  result  is 


JTM  -  u«  £ 

r+a  or 


(A. 6) 


with 


(A.7) 


R  -  (x2  +  x2)l/2 


(A.8) 


\2  a  +  2ra  cos9 


(A.9) 


To  evaluate  we  introduce  the  variable 


z  *  tan  9/2  ,  -  »  <  z  <  ® 


(A. 10) 


so  that 


q  l-z^  2dz 

cos  9  - - -  ,  d$  = - - 


(A. 11) 


It  is  readily  shown  that 


♦  s  ♦o  + 


(A. 12) 


where 


♦  =  p  (i-zz>-  <iz 

0  (ra) J_m  (l+z2)2  \a2+zz/ 


(A. 13) 
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t  ?  wi ,  >  v  j*t+*-*t  rf>  /?  ft** 


♦l  =  9(x)A  J 


_ (l-z2)(A-z2)dz _ _ 

(z2+1)^2(z2+B2)^2(AZ+z2) 


(A. 14) 


so  that 


g(x)  = 


(r-a)[x2+(r-a)2]1/2 


(A. 15) 


A(A2+1)  .  1  i 

A'1  (1«2)(Z2+A>J 


(A. 16) 


rx2  +  (r,a)211/2 

1x2  + 


(r-a)‘ 


(A. 17) 


^2)  ,lt  z  .  .t...  Al*?iU  •  _i 

.  *  1  »  i  A  A  1  r% 


(1+z2)(z2+A2)  A-1  z2+1  A_1  z2+A2 


(A. 18) 


*{» 


°"(r-a)J.«  1«2  ^  (1+Z2)2 


(A. 19 ) 
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All  of  the  integrals  in  Equation  (A.19)  may  be  found  in  Reference  8,  pp.  13 
through  17.  The  final  result  is 


2  *  _  * 
=  (r-a)(A+l)  "  r 


(A. 20) 


Thus 


Lim 

r+a 


n 

a2 


(A. 21) 


To  evaluate  we  first  note  that 


1  =  A-l 
r-a  2a 


(A. 22) 


Thus 


1  "  a[x2  +  (r-a)2]l/2 


G{A,B) 


(A. 23) 


where 


G(A,B)  =  (A-l) 


(z2+l)(z2+A2) 


_  dz _ _ 

(z2+32)l72(z2+1)l/2 


(A. 24) 


or  after  using  Equation  (A. 18)  we  can  write 

G(A,B)  =  (A-1)Gq  +  2GX  -  A(A2+1)G2  (A.25) 


where 


and 


(A. 33) 


where  K(k)  in  the  standard  elliptic  integral  of  the  first  kind;  i.e., 


K(k)  = 


'/l-k^sin^i 


(A. 34) 


In  the  complex  z  plane  the  integrand  of  G2  has  branch  points  at  z  =  ±  i 
and  z  =  ±  iB  and  simple  poles  at  z  =  ±  iA.  With  the  contour  shown  in  the 
sketch  we  can  replace  G2  by  an  integral  around  ;he  branch  cut  and  the  pole  at 


7  ~  1  A.  The  result  is 


_ _  rs 

A(a2-1)1/2(A2.b2)1/2  *  J, 


_ dt _ 

{A2-t2)(t2-l)l/2(82-t2)1/2 


(A. 35) 
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where  the  first  term  is  the  residue  at  z  =  iA.  The  integral  is  given  in 
Reference  8,  p.  41,  so  that 


G2  =  - 


it 


A(A2-1)1/2(A2-B2)1/2  BA2 


K(k) 


BA2(A2-1) 


n  - 


A2k2 

A2-l 


'(A. 36) 


where  n  is  the  elliptic  integral  of  the  third  kind;  i.e.. 


n(p,k) 


dt 


( 1+pt2  )\/(  l-t2 )  ( l-k2t2)  * 


0  <  k  <  1 


P  >  -  1  (A. 37) 


With  Equations  (A. 31),  (A. 33)  and  (A, 37)  we  can  write  G  in  the  following  form: 


G(A,B)  =  -  |  |l  +  jj  K(k)  +  K’(k)  (A. 38) 


The  final  step  in  the  derivation  of  the  ring-wing  kernel  is  to 
differentiate  Equation  (A. 33)  with  respect  to  r  and  take  che  limit  as  r  *  a. 
It  i«  useful  to  note  that 


Lim 

r+a 


AI 

dr  A 


_L 

2a 


(A. 39) 


and 


Lim 

r+a 


(A. 40) 
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Also 


and 


Lim  n 
r+a 


ak 


ar 


r=a 


x2/(x2+ 4a2)3/2 


where  E(k)  is  the  elliptic  integral  of  the  second  kind;  i.e.. 


(A. 41) 


(A. 42) 


E(k) 


n/  l-k^si  n24> 


d4> 


(A. 43) 


After  a  considerable  amount  of  detailed  algebraic  manipulation,  the  complete 
ring-wing  kernel  function  can  be  expressed  in  the  following  form: 


m.(i .  k\2 
1 


(A. 44) 


where 


k  = 


_ 2a _ 

V^+4a2 


0  <  k  <  1 


(A. 45) 


We  conclude  this  appendix  by  stating  two  asymptotic  properties  of  the 
kernel.  First,  we  note  tnat 
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and  the  kernel  has  the  simple  approximation 

K(x)  ~  H(x)  for  | x {  »  2a 
a2 


i 

(A.47)  j 


if 


t 


where  H(x)  is  the  heaviside  step  function.  We  refer  to  this  asymptotic  form 
as  the  slender  body  kernel.  The  second  result  is  obtained  for  x  «  2a  or 
k  ■*■!.  In  this  case  we  identify  the  Cauchy  singularity  in  the  kernel;  i.e.. 


K(x) 


_2_ 

ax 


for  x  «  2a 


(A. 48) 


